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Using the framework of wedge dynamics, we compute the effective transverse mass of a soft quark 
mode propagating in the expanding background of hard quarks and gluons created at the earliest 
time of the collision. We discover that the wedge dynamics does not require any external infrared or 
coUinear cut-off. The effective mass is produced mainly due to the forward quark-quark scattering 
mediated by the longitudinal (in sense of Gauss' law) magnetic fields. Contribution of the radiation 
field is parametrically suppressed. 
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O ■ I- INTRODUCTION 
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In the first paper of this cycle ( further quoted as paper [I]), we formulated a program that might result in a 
theory of ultrarelativistic nuclear collisions which is free from collinear problems and naturally establishes the infrared 
ps) , boundary for the space of "final" states at the very early stage of the collision (< Ifm). We have demonstrated that 
^ I even at very early times (much less than is required for any kinetic process to develop) , the collective interactions in a 
t> ' dense system provide the final states of the QCD evolution with finite dynamically generated masses that shield mass 
I singularities in the evolution equations]^. It was shown also that the null-plane dynamics are incapable of describing 
] local screening effects, because any type of kinetics is frozen on the light cone. It was suggested, that a more adequate 
[ approach requires the change of the global Hamiltonian dynamics which is used for the field-theory description of 



o 
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nuclear collisions. We proposed the so-called wedge dynamics which employs the proper time r measured from the first 
touch of the Lorenz-contracted nuclei as the Hamiltonian time. Our initial estimates in paper [I] were very qualitative. 
Clnj In two consecutive papers (further quoted as papers [II] and [III]), we have studied, in detail, the space of states 
O i' of wedge dynamics. In paper [II], we extended the qualitative analysis of scalar fields initiated in [I], and found that 
fH for the charged fields, the early-time evolution of the wave function is accompanied by a gradual rearrangement of 
^ the charge distribution, starting from its almost uniform spread along the light cone at r — + 0, and up to a narrow 
[ wave packet with a well defined rapidity at later times. We have shown that this re-distribution of the charge leads 
H " to currents in the rapidity direction and that these currents are the largest at the earliest r. The magnetic fields 
■ ■ ■ ' generated by these currents can be responsible for the interactions between the currents at the earliest moments of 
the QCD evolution. In paper [II], we studied the states of fermions in wedge dynamics and found the fermion field 
correlators that are used below for the calculation of the quark self-energy in the expanding system. In paper [III] , 
we addressed the issue of gauge fields in wedge dynamics. Several important problems were solved there. The natural 
gauge condition of wedge dynamics, — 0, was proved to be completely fixed (at the level of perturbation theory). 
The second (technically nearly most difficult) problem solved in paper [III], was the separation of the longitudinal 
(i.e., governed by Gauss' law) field and the field of radiation. In that paper, we also quantized the gauge field in the 
scope of wedge dynamics and explicitly found the Wightman functions and retarded propagator of the gluon field 
which are used in this paper for the practical calculation of the fermion self-energy. 



^ The idea that screening effects should be taken into account at the early kinetic stage of a collision has been articulated 
earlier and with different motivations by Shuryak and Xiong M and by Eskola, MuUer and X.-N. Wang M 
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Our decision to begin the exploration of potentialities of the wedge dynamics with the computation of quark self- 
energy is motivated only by technical reasons. The gluon propagator of wedge dynamics is a very complicated function, 
and we preferred to start with the computation of the fermion loop which has only one gluon correlator in it. We hope 
that the possibility of a technical simplification (compared to what we had to start with) discovered in this paper, 
will allow us to address a more important problem of the gluon self-energy in a reasonably economic way. 

In the course of this study, we employ a single heuristic assumption (supported by the analysis of paper [II]) that the 
field states with large transverse momentum, even at very early times, may be associated with the localized particles 
and thus can be described by the distribution with respect to the rapidity and transverse momentum. Our strategy 
of looking for the leading contributions, as well as all our approximations, in the calculation of the material part of 
the quark self-energy are based on this assumption. If it appears incorrect, then it is most likely that the quark-gluon 
matter created in the collision of two nuclei never and in no approximation can be considered as a system of nearly 
free and weakly interacting field states. 

II. FERMION RETARDED SELF-ENERGY 

In order to find the normal modes of the quark field in the expanding quark-gluon system , we are going to solve the 
Dirac equation with the radiative corrections, which can be derived as a projection of the Schwingcr-Dyson equation 
for the retarded quark propagator onto the one-particle initial state. For the quark field without Lagrangian mass, 
this equation reads 



i7^(xi)V^(a;i)-0(a;i) = J d X2^[ret]{xi,X2)ip{x2) ■ (2.1) 

The covariant derivative (x) of the spinor field in the curvilinear coordinates of the wedge dynamics includes the 
spin connection and it was found explicitly in paper [II]. For all calculations below, we employ the mixed representation 
which is the most profitable in heavy-ion problems. We are looking for the radiative corrections to the wave function 
with a given transverse momentum pt and rapidity 9 with the expectation that within the rapidity plateau nothing 
will depend on 9. However, we cannot totally eliminate the coordinate i] from the theory. We have to keep it 
explicitly, since the problem of the expanding field system cannot be reduced to (2+1) dimensions. In its expanded 



form, Eq. (2J) reads 



^OTi 2ti' Ti orji 



rTi poo 

/ dT2 / T2d?72S[ret](Pt;Ti,T2;?7i 
Jo J -oc 



V'(Pt;n,??i) = / dT2 T2dr]2^ret]iPt;n,T2;r]i~V2)-4'(Pt;T2,r]2) ■ (2.2) 



The retarded self-energy is an object that naturally emerges in the Schwinger-Dyson equation for the retarded prop- 
agator in Keldysh-Schwinger formalism Below, we employ its modified form developed earlier with the view of 
application to the inclusive and transient processes. We employ the notation used in Refs. [^HJl^.^ In this notation, 
the one-loop retarded fermion self-energy in coordinate form is 

■ 2 

2 

The two subprocesses that contribute this self-energy are depicted below. 



S[.et](xi,X2) = '-^rrGlret]ixi,X2)t'j'D\^^^^{x2,Xi)+t'^J^G[^ix,,X2)t'^^ . (2.3) 



^The indices of the field correlators with the Keldysh contour ordering of the field operators (like G[ab] ) as well as the labels 
of their linear combinations (like G[ret]) are placed in square brackets. 
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FIG. 1. The retarded forward scattering amplitude is contributed by two subprocesses, qg —> qg and qq —> qq. 

The retarded and advanced quark and gluon propagators G[ret] ^nd i'j^j,] were found in papers [II] and [III] of 
this cycle and are connected with the commutators G[o] and Z?[q] , 

G[,et](xi,X2)=0(Ti-T2)G[o](xi,X2), DIZ,]{x2, Xl) ^ -Bin - T2)DI^ {X2, Xl) + D'^' {x2, Xl) , (2.4) 

where D^l"{x2, xi) is the longitudinal part of the gluon propagator (governed by Gauss' law), and it enters in Eq. (2.4) 



in such a way that the condition D^^et] ^D^adv] = ^[o] is satisfied and the non-causal longitudinal part of the propagator 
does not violate the causal properties of the commutator Z?[o] • The correlators G[i] and D^^^ include densities of vacuum 
states as well as the information about the occupation numbers (phase-space population). Eventually, we shall prove 
that an approximation of the boost-invariance (infinite rapidity plateau) is not corrupted by any kind of cut-offs (the 
vacuum part never is). Therefore, all correlators (G, D, and E) will depend on two times ri and T2 separately, the 
difference of rapidities r] = rji — r]2, and the difference r = r 1 — r 2 of distances in xy plane. The latter is Fourier 
transformed to the transverse momentum dependence. In this mixed representation, 

^[ret]iTl,T2;r],pt) 



2(27r) 



d^fct[i"7™(n)G[„,](Ti, r2; 77,pl + fcOiV(T2)i^fil,™(r2, n; -ry, fct) 
-Ki"7™(ri)G[i] (ri, r2; 7y,pl + fcOtV(T2)i?['l„],™(r2, n; h)] , 



(2.5) 



where 7''(t) — j^/t. As it has been shown in [I], all fermion correlators G^^] can be decomposed as 
Gh (ti, T2; r/, qt)^qt [.g[L]7° + Sht'] + g^^rY + i5[l]g.e™7"75 

where, for the sake of brevity, we denote eft = pi + kt. A similar decomposition takes place for the self-energy, 

^ret]{TuT2;v,Pt) = + S^t' + ^'^ QrY + iSVe™7"7' 



(2.6) 



and we obviously have. 



i(±) 
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It becomes easier to analyze the various pieces of the quark self-energy if the gluon correlators 
the form of the following decomposition]^ 



(2.7) 



(2.8) 



are taken in 



^In what follows, we use the Greek indices for the four-dimensional vectors and tensors in the curvilinear coordinates (index 
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('J'?) 



n, , - ^-r,('-')) 



Dr 1 - tlT)^^'') 



(2.9) 



where the first term in D[a]rs is due to the transverse electric mode, and all invariants of V^adv] (except T^^'^j^^ ) have 



two terms, I?L ^ from the transverse magnetic mode of the radiation field, and 2?^ , from the longitudinal field. All 
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(...) 



ong\ 



these components were found in paper [III] and are given in the Appendix in the form which is used in the calculation 
below. After some algebra, we can present the retarded quark self-energy in the form. 



^[ret]iTl,T2;r],pt) 



2n 



where the scalar invariants of S[r-et] are the bilinears of the fermion and gluon scalars. 
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(2.10) 



(2.11) 
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T(±),^(TE) ^(2). 



jqtPt) T(T)^(m) 
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Tl 



T2 



In these equations, the sum X][q b] runs over [a, f3] = {[ret, 1], [1, adv]}. 



III. FERMION MODES IN THE EXPANDING SYSTEM 



(2.12) 



We shall look for the dispersion law of the fcrmions in the proper-time dynamics studying the Dirac equation (2.2) 
with radiative corrections. Since the fermions are massless, it is convenient to use the spinor basis where the Dirac 
matrices are 

7° 

and the Dirac equation can be split into two separate equations for the left- and right-handed two-component spinors. 
The latter reads as 

1-1 , 
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T2d-q2^R{puTi,T2]in ~ m)i^R{puT2,'n2) , (3-1) 
where the matrices of the right-handed differential operator G]^^ and of the right-handed self-energy E/? are 

r ;/a , 1 

Gr^{Put, rf) 
^R{pt;Ti,T2;r]i - 772) = 



«(5r + ^ - ^djj) p^ - ipy 

Px + iPy i{dr + ^ + ^d.,j) 

S^(') -(p.-zp,)S^(+) 
-(p,+ipj,)E^(-) 



(3.2) 



77 is an exception, it always denotes the rapidity direction), and the Latin indices from a to d for the vectors in flat Minkowski 
coordinates. We use Latin indices from r to 10 for the transverse x- and y-components (r, ...,w — 1, 2), and the arrows over the 
letters to denote the two-dimensional vectors, e.g., k = (k^, ky), \k\ = kt. The Latin indices from i to n (i, n = 1, 2, 3) will 
be used for the three-dimensional internal coordinates — {x,y,r]) on the hyper-surface r = const. 
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The equation for the left-handed spinors differs from ( |3.l| ) only by a change of some signs in matrices ( |3.2| ) and leads 
to the same dispersion law. A solution with positive energy is looked for in the form, 

i^niPuO;r,,) = ( J^^^^^ ) e--™-(-^) , (3.3) 

where fi is the effective "transverse mass" of the mode. For the free on-mass-shell solution we have ^ = pt- To solve 



Eq. (3.1), we introduce an auxiliary (left-handed) spinor 

(-(-^')/.) ■ (3-4) 

We insert ( |3.3| ) into ( |3.l| ), multiply it from the left by (^^) and integrate this along the hypersurface ti — const. 
Then the left side of the equation becomes 

Tidrii4!{pt,6';Ti,7ji)G^^{pt;Ti,rii)'iljR{pt,9;Ti,f]i) = An- — — pi 5{e - 9') . (3.5) 

In deriving this equation, we assumed that /i is independent of ti. The weak dependence is admissible, provided 
dfx/dri <^ a solution that has this property does indeed exist. The right hand side of the equation is, in fact, 

independent of 9' and is of the following form. 



PTI POO 

pI / dT2 / TiT2dr]2d9d{r]i - 772)6 

Jo J -oc 



ifiTi cosh(r7i +^) z/ir2 cosh 772 



X [e 5 — + e — 5 — + e ^ — 1]^+ + g — ^ — 1]^-] , (3.6) 

where the exponentials are due to the Thomas precession of the spinor field. Next, we integrate both sides with 
respect to 9. Two rapidity integrals, d9dri2, on the right absorb the precession factors yielding the product of Hankel 
functions, 



7T'H[%i^T,) H[%{pr2) = -^=e^^^^^-^^^ . (3.7) 



Finally, we arrive at the dispersion equation that defines the fermion "transverse mass" /i as a function of the transverse 
momentum and the latest time ti. 



^ Jo J -00 



(3.8) 



As it has been discussed in paper [II] for fermions (similar arguments are true for gluons), only the independence 
of the quark and gluon occupation numbers nj and ng on rapidity can provide that the invariants S*^*-^^ and S"^^^"^ 
naturally depend only on the difference rj — j]i — ri2- We shall consider only this case of the local homogeneity; we can 
do it safely only because no coUinear singularities which may require a rapidity cut-off (e.g., ±1") in the phase space 
will appear in the theory. Since we are computing an essentially local quantity, such a cut-off would be unphysical. 

as 



With this reservation, we may rewrite Eq. ( 3 



M(Pt) ^Pt+ dT2V^e'''(f')(^^-^^)[E"(ri,r2) +ptE^(Ti,T2)] , (3.9) 







where we introduced the notation, 



E"(n,r2) = ^ E / ^^fc* r ,.9UPr,,^'^< + T;^^??]' (3.10) 
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S^(ri,T2) 



[",/3] 



- 2 



,^(TB) „(2). _(9tPt)_^(w) 



[/3] 



(3.11) 



Comparing these equations with Eqs. ( ^.ll| ) and (2.12), we may observe a significant simplification. The terms with 
the off-diagonal components and T)^"!") have dropped out. These terms, as it can be seen from Eqs. ([A^)- 



(A7), (AlO)-(All), and ( Al§| )-( pVl9 ), are odd with respect to r/, while the invariants 5° = g 



L{+) 



1'^^ ) and 



are even. Therefore, integration over rj eliminates the terms with the off-diagonal components. 



IV. PROPAGATORS, DENSITIES OF STATES, AND OCCUPATION NUMBERS IN THE EXPANDING 

SYSTEM 



In this section, we collect condensed information about various correlators of quark and gluon fields derived in 
papers [II] and [III] which are necessary for the calculation of the quark self-energy. We also discuss our specific 
choice of occupation numbers ng{kt,a) and nf{qt,9). All field correlators are defined as the expectation values over 
the distribution of the background particles. The latter are the excitations of the modes allowed by the constraints 
and the boundary conditions of wedge dynamics. The Fock space of these excitations was constructed in papers [II] 
and [III]. We have analyzed two sets of quantum numbers that may label the states. Both sets include the transverse 
momentum pt and polarization index. In one set, the remaining variable was the boost v (the variable conjugated to 
the coordinate 77); this set proved to be very useful in the practical calculation of the gluon propagators. In the second 
set, the particles are labeled by their velocity Vz = tanh^^ in the direction of the collision axis. This representation is 
used below. The fermion spectral functions are 



G[io]{qt,0;Ti,T2) = [1 - n+{qt,9)]G''^l^{qt,0;Ti,T2) - (g*, 6I)G[qj'] (gt, 6*; ti, T2) , 



(0) 



G[oi](gt,6';Ti,r2) = (gt,6')G'[°o](gt, 

Their gluon counterparts are of a similar form. 



;n,T2) 



[l-n7(<7t,0)]Gg](gt,0;ri,T2) 



(4.1) 



D[iQ]^{kt,a;Ti,T2) = [1 + ng(fct, q;)]!)?^ (fct, a; n, r2) 



where Z3^°j and Gj^J are the vacuum correlators of a given type [a]. They are defined as vacuum expectation values 
of the binary products of field operators. 



ng(gt,Q;)Z)[°^](gt,a;Ti,T2) , 
D[oi\{kt,a;Ti,T2) = ng(fct, a)!?!"^^ (fct, a; n, r2) + [1 + ng(gt, a)]D[p^] (g*, a; n, T2) , 
.(0) 



(4.2) 



Gf?l,(xi,X2) 



m^{xi)^{x2m , G[°j](xi,X2) = i{Q\^{x2)^{xi)\0) 
Xxi,X2) = -i{Q\A^{x2)Ai{xi)\Q) 



D 



(0) 



{xi,X2) 



'[10]/mV-^l>-^2; - -i{Q\Ail{xi)A{x2),nP) , 
In this approximation, the field (anti-)commutators. 



D 



(0) 
[01]/? 



(4.3) 



^[0] — G[io] — G[oi] — G^ni — Gfnii — G 



Du 



'[0] - -D[io] - -^-'[oi] - --[10] - ^[01] - -^[0] 
appear to be insensitive to the presence of the particle distribution, while their counterparts, 

- Gno] + Groii = [1 - -1^^°^ 



Du 



D 



(0) 

[10] ^ 

(0) 



D 



(0) 
[01] ■ 

(0) 



D 



(0) 

[0] ' 

(0) 



(4.4) 



= i?[io] + i?[oi] = [l + 2ng]D'^^ 



(4.5) 
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include the occupation numbers which modify the original vacuum density of states. For the sake of simplicity, we 
take riy — nj = Ji/, which corresponds to a neutral system. 

The Wightman functions ( |4.l| ) and (4.2) (or their various linear combinations G[^] and eventually appear 



under the integrals dO and da. One must keep in mind that in order to reduce and oi?} to the standard form of 



the vacuum correlators, at least two shifts of the integration variables is necessary. Only after that will Gj°j and -Dj^j 
explicitly depend on the boost-invariant variables i] and T12. The functions ng{kt, a) and nf{qt, 0) are not indifferent 
to this shift. It may well happen that a formal shift in 6 or a will drive the stationary points of the wave functions 
or the singularities of the field correlators outside the physical boundaries of the distributions ng{kt, a) and nf{qt,6). 
Therefore, different representations of G[i] and -D[i] must be used for the study of different subprocesses. One has to 
account for the reservations stemming from the derivation procedure described in Sec. 4 of paper [II]. These different 
representations of the quark and gluon correlators are quoted in the Appendix. 

In our picture, first outlined in paper [I], the fermion vacuum mode with small transverse momentum pt and zero 
rapidity is modified by its forward scattering either on gluons with high momentum kt and rapidity a, kt 3> pt or on 
quarks with high momentum qt and rapidity 6, qt ^ pt- These hard modes are created at the earliest moment of the 
collision and can be treated as well formed particles by the time t ~ since at that time rkt 3> 1, and rqt 3> 1. 

Therefore, they may be consistently described by the distributions, 

f m 0iQt~P*) . J^g 0{kt~p.) 

nfidue)^^^^, n,(fc„a)«-^^^. (4.6) 

where is the lower bound of the "hard" partons distribution. Both distributions (per unit area, per unit rapidity) are 
chosen on purely dimensional grounds, since we believe that the creation of a parton with large transverse momentum 
is described by perturbative QCD which has no intrinsic scale. 

Currently, the normalization factors Afg and A// are the only (apart from the coupling as) parameters of the theory. 
The cross section ttR^ and the full width 2Y of the rapidity plateau are defined by the geometry of a particular 
collision and the c.m.s. energy, respectively. These are irrelevant for the local screening parameters we are interested 
in. In the first approximation, one may try to extract them from the event-by-event measurement of the high-pt tail of 
the collision products and incorporating the standard phenomenology of the fragmentation functions for the analysis. 

As it was pointed out in paper [I], even in dense systems, the QCD evolution at large is not likely to be 
affected by finite-density effects. Thus, one may also try to employ the known structure functions (without shadowing 
corrections) and the factorization scheme in order to estimate Mg and A//. A most appealing opportunity to find 
ng{kt,a) and nf{qt,0) from first principles, associating them with the known properties of hadrons and the QCD 
vacuum, is still very distant. 



The distributions (4.6) are used below with the following informal reservations. First, the total energy of any 
collision is finite and kt and qt have (though very high, but finite) upper boundary. Eventually, this leads to the 
self-energy which is free from collinear singularities in the interaction of charges with the vector gauge field. Second, 
though the distributions (4.6) are boost-invariant, only the particles which physically affect the forward scattering 
must be accounted for. There is a strong correlation between the position 77 where the particle with large transverse 
momentum qt is measured (or is interacting) and its rapidity 9. Hence, the limits of integrals da and dO over the 
rapidities of real particles (which either mediate the scattering or are in the final states) cannot exceed the actual 
rapidity boundaries of the scattering process. In its turn, this puts an additional requirement on the notion of the 
distribution itself. It must be normalizeable in the physical volume of the reaction. This volume is defined, in fact, 
by the light cone (i.e. causality of the forward scattering amplitude). [We remind the reader that the notion of a 
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distribution itself makes sense only after it is prepared (measured) at least in a gedanken experiment. Hence, the 
distributions ng and nf must exist, in this sense, both at final time ti and at the initial time T2 in the expression for 
the self-energy. In its turn, this limits the time T2 from below.] 

V. LEADING PART OF THE DISPERSION EQUATION 
A. Derivation of the dispersion equation 

The most important outcome of this work is that the major contribution to the effective quark mass comes from 
the ?7r7-component of the propagator of the longitudinal field. This contribution is computed in all details below. All 
other terms are associated with the propagation of the transverse fields and they appear to be parametrically small 
in the domain Tipt < 1, (ti — r2)pf <C 1, ti — T2 <C ti, where the dynamical mass of the fermion is effectively formed. 
The component I?,,,, of the propagator establishes the connection between the A,-/ component of the potential and 
the jri component of the current. In its turn, A,, is responsible for the ry-component i?^ — drA^, of the electric field 
and the x- and y-components, = dyAjj, By — —dxAjf of the magnetic field. The electrical field in the longitudinal 
77-direction is not capable of producing scattering with transverse momentum transfer. However, this transfer can be 
provided by the magnetic forces; the two currents can interact via the magnetic field Bt ~ [Bx, By). The origin of 
these currents is intrinsically connected with the geometry of states in the wedge form of dynamics. Any state with 
a given pt begins its life being widely spread along the light cone. If the state is charged, then local charge density is 
small. With time going on, the spread of the wave function diminishes and the charge become localized in a narrower 
rapidity interval (see Fig. ^). Therefore, any charged state carries a current in the longitudinal (rapidity) direction. 




FIG. 2. Evolution of the charge density in the typical state of wedge dynamics 

The magnetic fields of the transition currents provide scattering with the most effective transfer of the transverse 
momentum. Indeed, at time T2 a quark with the transverse momentum pt, T2Pt <SC 1, interacts with the gluon field 
and acquires a large transverse momentum kt, T2kt 3> 1. This transition is characterized by a drastic narrowing 
of the charge spread in the rapidity direction, and must be accompanied by a strong 77-component of the transition 
current. A similar transition in the opposite direction happens at time ti, when the gluon field interacts with another 
quark that has large initial transverse momentum kt, and recovers the soft state with T2Pt ^ 1 in the course of this 
interaction. This second transition current readily interacts with the magnetic component of the gluon field. Our 
estimates indicate that the leading contribution comes from the term of D^^(t2, ti; 772 — ?7i; fct) which is proportional to 
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5(r]i —772) and does not depend on kt [in coordinate representation, this term is just proportional to S{r]i — 772)^(^1 ^1^) 
]. This is a long-range contact interaction of the two currents, and is not limited by the light-cone boundaries (which 
suppress the interaction via the strongly localized states of the radiation field). Furthermore, the contact part of 
the longitudinal propagator is the only one that brings into the integrand of the dispersion equation ( |3.9D the term 
which is singular at ti — T2 — s- 0. Therefore, it is capable of providing an appreciable contribution into the effective 
quark mass, which is defined locally. This part of the self-energy allows for an exact calculation with a simple analytic 



answer which is presented below. Estimates of all small terms are explained in Sec. VIl and VIIl . This contact part 
of the longitudinal propagator is 

4~'(^2,Ti;r72 - Vi;kt) = ~^^^Si^). (5.1) 
Because of the extreme locality of provided by S{ri), the invariants g^i^ of the fermion density function in 



(3.10) lose their kinematic coefficients. 



First, we integrate over r/, which leads to = ('''1 ~ ■''2)^- Next, we change d'^kt for d^cft and integrate over the 
orientation of ql gaining the factor 27r in . In S^, g'^-^ is integrated with the weight factor {qt ■pt)/pf. Therefore, this 
term identically vanishes after integration over the azimuthal angle. The only remaining integral over the transverse 
momenta of hard partons is 

Yi[[Ti - T2)qt\ dqt = . (5.3) 

Tl - T2 



P. 



Eventually, we may write the dispersion equation (|3.9|) as follows 



M - + r d.2^e-(— )yo[(n - rM . (5.4) 



B. Study of the dispersion equation. 

According to the qualitative analysis of paper [II], the dynamics of states is different in the two limiting cases, 
Tpt < 1 and Tpt > 1. With respect to pt ^ l/r, the states are divided into "hard" and "soft" states. Therefore, it is 



natural to take ~ pt in Eq. (5.4). Further, it is convenient to trade variable T2 for y — (ti ~ T2)/ti, 



Kp^ . 1 + ^-^^frm r ,yi_iae-(--o-yo(np.,) . (5.5) 
Pt ^T^RtPt Jo V^-y 

In this form, the dispersion equation clearly reveals two distinctive regimes. When rpt < 1, then the function Yq{x) 



behaves as a logarithm, and the right hand side of Eq. (5.5) becomes proportional to \n{2 / Tipt) , the effective width 



of the rapidity interval occupied by the state at the early time of the evolution. When rpt > 1, then Yo{x) ^ l/^/x, 
and the integral becomes proportional to l/^rpt, the effective rapidity width at later times. Thus, the dispersion 



* The authors appreciate discussions with Edward Shuryak, who pointed out that the small effect of the radiation fields is 
much less surprising than the finite contribution from the longitudinal fields. 
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equation ( ^.5| ) clearly reveals two distinctive regimes which were qualitatively analyzed in paper [II] . The solution of 
Eq. (5.5) is generally complex. Taking /i = /i' + i/i" we can separate real and imaginary parts of this equation, 



Ti/i - TiPt 



f 



dy 



1 



dy 



y 



sm{fj.'Tiy)Yo{Tipty) 



cos{n'Tiy)Yo{Tipty) 



(5.6) 



(5.7) 



(The unit upper limit in these integrals corresponds to T2 — 0, and is, as a matter of fact, fictitious. Practically, we 
are interested only in the domain where T2Pt ^1.) We have rearranged the factor in front of the integral in such a 
way, that at early times, this factor is small. It has been shown in paper [III], that the longitudinal part of the gluon 
propagator vanishes when the distance rt exceeds ri. Therefore, this factor is proportional to the (small) number of 
hard partons per transverse area occupied by the soft quark mode. Hence, we can analyze Eqs, ( ^.6[ ) and (5.7) by 
successive approximations. It is clear, that in the lowest approximation, we can take /i' — pt in the RHS of these 
equations, and that the imaginary part p" can be neglected. Using 



cosa;~l, Yq{x) ^ 271 "^[7£; + lna;] 



as an approximation, and computing the remaining integrals, we arrive at 



Pt T^^Rfpf 



Pt 



TlPt 



IE 



In- 



TiPt) 



In- 



TlPt 
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(5.8) 
(5.9) 



These dependences are plotted below as functions of the argument Tipt up to the pre- factor asCpJ^f /tt^R^pI 




FIG. 3. Self-energy corrections to the real part {^i' /pt — 1, upper curve) and to the imaginary part [p" /pt, lower curve) of 
the effective transverse mass as functions of npt. 



Eq. (5.5) describes the evolution of the effective transverse mass /i of the state with a given transverse momentum 
Pt as a function of the proper time ri. We see, that the real part /i' gradually grows with time reaching its maximum 
at TlPt ~ 1- The mode acquires an "adjoint mass" due to the interaction with hard partons, as was anticipated. The 
curves cannot be trusted above the boundary ri « since at later times, the mode becomes "hard". It cannot be 
viewed as a soft cloud swept with uniformly distributed hard particles. The condition 

dri T\ ' 
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which was used in the course of the dispersion equation derivation, is clearly fulfilled near the maximum of the 
dispersion curve. 

One more important dependence is hidden in the pre- factor agCp.N'f /i^^R'tPti and is not visible from the figure 
above. This factor scales as p^^, clearly indicating that at large pt, the effect of screening is small. There is almost 
no hard particles with kt.qt > Pt- 

VI. CANCELATION OF COLLINEAR TERMS IN THE VACUUM PART OF THE SELF-ENERGY 

Usually, the self-energy E is studied in the momentum representation, and the first subject of concern is the 
ultraviolet divergence of this function. It is well known that this divergence can be at most logarithmic. Thus, when 
we compute the four-dimensional integral over the momentum in the loop, this divergence can show up only in the 
last of these integrations. We compute the self-energy in the mixed representation. Hence, we cannot see the UV 
divergence explicitly but we must already see (e.g., in ImE) the various infrared divergences which emerge due to real 
processes with massless fields. A corresponding analysis for the case of the null-plane dynamics was done in paper 
[I]. These divergences must be regularized (or even removed from the theory, as is done by means of dimensional 
regularization) before the UV renormalization. The primary goal of this section is to demonstrate that in the wedge 
form of dynamics the quark self-energy is completely protected from coUinear problems, and that this is not a surprise. 
Indeed, in the theory with massless fermions and gauge bosons, the infrared singularities show up in a different way 
depending on the type of Hamiltonian dynamics (including the gauge condition) which is used to describe the process. 
In the gauge — they look like coUinear divergences. In the gauge = 0, they look like an infrared problem of 
the proper field of the charged particle. In both cases, the problem emerges due to the incomplete gauge fixing, and 
manifests itself through spurious poles of the gauge field propagators. As has been shown in paper [III], the gauge 
A'^ — is fixed completely, and therefore, the quark self-energy that we compute here is totally free of these problems. 

In order to demonstrate this appealing feature we shall compute (the most dangerous in this respect) the vacuum 
part of the fermion self-energy, concentrating on the terms where the integrand as a function of the rapidity a is not 
suppressed at |a| ^ oo . A self-consistent piece of this type is the contribution of the transverse electric mode of the 
gluon field. The tensor part of any gluon correlator for this mode is of a very simple form, 

D^,,^i6rs-kr.h/k^)Vlf ; (6.1) 

it has no rj components, and the scalar functions T^fJ^ can be computed exactly since they have simple integral 
representations. We use this piece of the self-energy to explain the principles we base our calculations on. Using 



Eqs. (|3.10D and (|3.1lD , we get 



[4i?]rre(^i'-2;ry,ft) = ^0(r^,)0(n-r2) J d%[q,g^^^\q,)Vf,f (h) - q,g;i^\q,)Vf,f (h)] , (6.2) 



where eft — kt + pt and the minus sign in the second term is due to the definition (2^) of D^adv]- The E^^^j is fully 
confined within the light wedge t|^2 > 0- Then, the vacuum quark and gluon correlators have the following form, 

I'fSf (^2, n; ry2 - ryi; h) = 2-^0{Ti2)Jo{Ti2kt) , 2?^f (t2 , n ; r/2 - m; h) = 2-^i Yo{Ti2kt) , (6.3) 



L(±) _ . Tie^" ~T2e " ^,_2 , , , ^ _ /, 2 L(±) _ T^e^" - r^e 



[0] ~ A^/YTI 



Oirt^) Mqt,J\Th\) , = , /rfi ^ifeVl^i^l) , (6.4) 



4^ "^--"^-'v— 4^ 
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where = '''i + '''I ~ 2riT2 cosh?], and the commutators Dj^^ and G[o] include the causal 6'(r^2) by definition, and 
the terms proportional to ^(— t^2) the densities Dj^^ and G[i] are omitted. Integration over the angle between 
the vectors pt and kt involves only the invariants g^^y According to Eqs.( |6^ and ( |6.4[ ), we have to integrate, 

qtJi{Ti2qt)dip = 2ti [kfh{Ti2Pt)Ji{Ti2kt) + PtJi{Ti2Pt)Jo{'ri2kt)] , 



/■2-7r 

Jo ' 



qfyi{Ti2qt)d(p = 2iT {0{kt - Pt)[ktJo{Ti2Pt)Yi{Ti2kt) + Pt Ji(Ti2Pt)Fo(n2fct)] + 
+9{pt - kt)[ktJi{Ti2kt)Ya{Ti2Pt) + PtYi{Ti2Pt)Jo{Ti2kt)]} 



(6.5) 



This integration is done with the aid of the so called addition theorems g] for Bessel functions of the argument 
qt — [k^ + + 2fcfPf cos^s]^/^. Starting from this point, we can continue in two ways. The most straightforward 
option is to use the gluon correlators in the integrated form of the Bessel functions, thus sweeping under the rug the 
singular behavior stemming from a oo. This leads to the integral 



(rhmri - T2) 



Tl2 



Ja{Ti2Pt) / k1[Ji{Ti2kt)Yo(Ti2kt) + J(i{Ti2kt)Yi{Ti2kt)]dkt -2ptJi{Ti2Pt) / ktJQ{Ti2kt)Yo{Ti2kt)dkt 



Pt 



+PtJi{ri2Pt) / ktJoiTi2kt)Yo{Ti2kt)dkt + Jo{Ti2Pt) / h JoiTi2kt)Yi{Ti2kt)dkt 



PtYl{Ti2Pt) / ktJoiTl2kt)JoiTl2kt)dkt -Yo{Ti2Pt) / fcj Jo(Ti2fcf)Ji(Ti2fci)dfct 



Pt 



(6.6) 



Here, all the integrals can be computed explicitly as indefinite integrals. The integrals from to pt yield the regular 
part of the answer below. Taking the upper limit of improper integrals to be A ^ oo, we get the singular part as the 
limit, 

A2 - 



lim 

A^oo T 



[jQiTl2Pt) + TptJl{Tl2Pt)]JliTA)Yl{TA) + Tpt Jl (ri2Pt ) -^0 (TA)ro (tA) 



(6.7) 



Using the asymptotic expansion of Bessel functions, we find that the singular part is built from S{ti2) and its derivative. 
The full answer reads as 

ISr„l./l = — — 6'(Ti2)6'(-ri - T2) 



167r 



T12 



^.2J2{TV2pt) ^ 

X i Pt ^ 

T12 



T12 



T12 



-^(ri2) 



4ptJl(Ti2Pt) 



T12 



S{ti2) 



(6.8) 



Thus, the answer is singular at the null planes rj^2 ^ 0. In order to find the true source of this singularity, we shall 
proceed in a different manner, keeping the gluon invariants 5^'-^^ in the integral form. We shall integrate over kt first 
and leave the integral over the gluon rapidity a for the end of calculation. This leads to 

J2(Tl2Pt) 



da 



Jo(n2Pt) 

+PtJl{Tl2Pt) 



ktJi{Ti2kt) cos{Ti2kt)dkt 



ktJo{Ti2kt) cos{Ti2kt)dkt 



Tl2 

CO 



27r Tl2 



ktYilTi2kt) s\n{Ti2kt)dki 



ktYQ{Ti2kt) sm(Ti2kt)dkt 



(6.9) 



where, we remind the reader that, T12 = Ti^T2 — T12 cosh(a — ?/;). The integrals in this expression are the well-known 
Fourier transforms of the Bessel functions, 
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kfJi{Ti2kt) cos{Ti2kt)dkt 



ktJo{Ti2kt) cos{Ti2kt)dkt 



kjYi{Ti2kt) sm{Ti2kt)dkt = 3Ti2Ti2 [T^2 " T\ 

/>OC 

/ ktYo{Ti2kt) sm{Ti2kt)dkt = -T12 [T12 
Jo 



-5/2 



2 1 -3/2 
'^12] + 



(6.10) 



where the distribution is defined in a standard way with the due number of subtracted terms of the Taylor 
expansion in the integral J f{x)x^dx ||l0| . Three different issues are important here. First, each of the integrals 
( |6.10 ) is a well defined distribution that includes all necessary regulators which provide the convergence of subsequent 
integrations. The Besscl functions themselves are defined as the Fourier transforms of the (+)-distributions and we 
just recover the original regular form by doing the inverse Fourier transform ( 6.10| ) (see Ref. |10|) . Second, as it 
will be shown in the next section, the singular behavior of the integrals ( |6.10|) originates from the collinear domain. 
The (+)-prescription that emerges here eliminates them term-by-term. Third, after the result of the term-by-term 



integration is put back into Eq. (6^), the singular collinear terms just cancel everywhere, including the null-plane 



T^2 = 0. This type of cancelation of coUinearly singular terms takes place in all other pieces of the vacuum part of 
the quark self-energy. 

All these observations lead us to the conclusion, that even in its vacuum part, the self-energy does not suffer from 
collinear problems, which seems to be a unique property of the expanding system. We do not continue to study the 
vacuum part of the self-energy here, since we are currently interested only in its material part which is discussed in the 
next section. [ The full analysis of this part, including the issue of its renormalization, will be published elsewhere. ] 

VII. RADIATION FIELDS IN THE MATERIAL PART OF THE SELF-ENERGY 



We found that the major contribution to the one-loop effective mass of a "soft" quark mode in the background of 
"hard" quarks and giuons comes from the quark-quark forward scattering mediated by the magnetic component of the 
longitudinal field. The purpose of this section is to demonstrate that the interactions via transverse fields (including 
the forward scattering of soft quark on hard giuons) is a secondary effect at least at the very early stage of the nuclear 
collision. This conclusion may look counter-intuitive, since, namely in the interactions of the transverse fields, we 



expect to encounter the collinear enhancement of the radiation amplitude. As it has been shown in Sec. VI, in the 
vacuum part of the self-energy, the integrals of this type (taken in the limits from to 00) cancel each other leaving 
the vacuum sector free from collinear divergences. The statistical weights Afg and Aff, which are different for the 
different terms, prevent such a cancelation in the material part. Thus we have to analyze each term of the material 
part separately. 



As in Sec. VI, we consider an isolated piece which corresponds to TE-gluons. The gluon correlators of this piece 
arc the most singular and are known not only in the integral representation, but in closed analytic form also. The 
last circumstance is very helpful for the analysis of the multiple integrals we meet below. (The terms identical to 
those computed below, also appear in the part of self-energy due to the TM-gluons; the remaining terms of TM-sector 
are less singular and, eventually, smaller than considered here.) The corresponding fragment of the quark self-energy 
( |3.10 ) in the dispersion equation (3.9) is 

liTE) 



S°(Ti,r2) 



[ri - T2 



dv <lt[9[op[i]^^ 



T,iTE)i 

f[i]^[o] J 



(7.1) 



S^(ri,T2) 



(TE) 



47r 



(ri - T2) 



d^h 



drj 



(pS) AktPt){ktqt) 



Pi 



ktpt 



(TE) 
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{TE) 

[0] 



(7.2) 
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The invariants of (anti-)conimutators are the same as in the vacuum case, 

2)cOsh(?7/2) 2^ J / N T r 

6I(ti2) Ji(gtTi2) , 5[o](n,T2;?7;(,,, ^ 



0/ ■ (n -r2)cosh(?7/2) 2 N 11 \ Ti cosh(ry/2) 2 N ji ^ ,^ o\ 



and the invariant 'D'EF{t2, ti; rj; kt) is given by the first of Eqs. (|6.3[). They all differ from zero only for the timelike 



[0] 

T12. Hence, the material part of the invariants and 2?^]^^'' will be needed only in this domain. As it has been 
discussed earlier, the distributions include only "hard" particles which are defined with respect to the soft mode with 
the transverse momentum pt by the inequalities, kt > p*, and qt > p*, where > pt- Now, we are interested only in 
the material part with occupation numbers given by the equations, 

and we must keep in mind the width 2Y of the rapidity plateau with the goal to study if this is a significant parameter 
for the calculation of local quantities. We may also question the validity of these formulae at sufficiently large kt and 
qt, since without a cutoff, the integral J dkt/kt diverges. 

The material part of the densities will be employed in two different forms, 

n , de , ^ , , (ri - r2)cosh(?7/2) , , , , , , 

g\i][Ti,T2;ri;qt) ^ / — n/(6'; gt) cosh6'smqtri2(6') = nf(qt)Yi(qtTi2), (7.5) 

^ ' J-00 T12 

rj f°° dd 

9[\\{Ti,T2]-q;qt) = - i cosh - / — nf{9;qt) cosqtTi2{d) = - icosh{rj/2) nf{qt) ^0(9*^12) (7.6) 
T^li]^\T2,Ti;r]]kt) ^ {niy^ j dang{a; h) cos ktTi2{a) = i Yo{Ti2kt) ng{kt) , (7.7) 



where the second equation in (7.5)-(7.7) is valid only when Ug and Uf are rapidity-independent, and we employ the 
following notation, 

Ti2{ct) = Ti cosh(a — 77/2) — T2 cosh(a + 77/2) — T12 cosh(a — ip) , 
T^2 = '''1 + '''2 ^ 2tiT2 cosh 77 > 0, tanhip{ri) — tanh — , 

Tl - T2 2 

I77I < 770 = In ^ « ^i— ^ = e, tanhi/^(±?7o) = ±1, 7/;(±77o) = ±00 , (7.8) 

T2 ^TlT2 

where ^ ~ (ti — T2)/ ^/TlT2 ~ 770 is the main parameter of our calculations. This parameter is supposed to be small in 
order that the notion of the current transverse mass /i(pt, ti) has the expected meaning of a slowly varying parameter. 
The geometric mean time t,„ = y^rprj has a simple interpretation. The two characteristics, one connecting the points 
(t2,— 770) and (ti,77o), and the second one connecting the points (t2,7;o) and (ti,— 770), intersect at the point (r„i,0). 
The proper time r„i is always inside the domain of the "causal interaction." 

Let us start the analysis of the radiation-dominated terms with the invariant [2*^]^^^ (The invariant [S"^]^j^^ 



appears to have an extra small factor ^.) According to Eqs. (7^), ( |7.5| ), and ( |7.7[ ) it can be written as a multiple 
integral. 



S°(n,T2) 



mat 87r2 



dvOir!2, 



/ dqt / dip Jo{Ti2kt) / d9 cosh 9 smTi2id) qt 



K (ti - T2) cosh 77/2 f°° dkt f^"" , , , [°° , , ^n, X , ^ 

—iST— — / / d^1tJi(Ti2qt) / da cosT^h U n - r2 , 7.9 

T^Rl T12 Jp, h Jo J_oo J 
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where we choose the integral form of the densities g[i] and f [j^^' in order to find the domain in the multidimensional 



space where the dominant contribution comes from. Since the two terms in (7^) are not expected to interfere (or UV- 
diverge), we are free to change variables in these terms independently. We leave (Pk^ in the second term, and change 
for (Pql in the first one. The next step is to integrate over the azimuthal angle between ql and pt in the first term 



of Eq. (6.2), and over the angle between kt and pt in the second term. This integration deals only with the retarded 
propagators inside the self-energy loop and selects the lowest angular harmonics, 

/•2Tr 

/ qtJi{Ti2qt)d(p = 2Tr [ktJo{Ti2Pt)Ji{Ti2kt)+ptJi{Ti2Pt)Jo{Ti2kt)] , h > pt , (7-10) 
Jo 

Jo{Ti2kt)dip = 2tt Jo{Ti2Pt)Jo{ri2qt) , qt > Pt ■ (7-11) 



Only the first of the two terms in Eq. (7.10), corresponding to the collinear geometry in the transverse plane survives 



in the limit of kt » pt and has to be retained by our major assumption. The second term describes the deviation from 



collinearity and is small. However, it is instructive to keep it for a while. After these angular integrations, Eq. (7.9) 
becomes 

-1 (TB) . ^ ^+oo r «/- 



S°(ri,T2) 



mat ^" 



/+00 ^ I'CO /■oo 

d'V^i'^i2)\—^ / dqtJo{Ti2Pt)JoiTi2qt) / d9 cosh 9 sin T12 (6) qt 
-00 Jp, J -00 

^^(n^Ts^cosh^ r ^[ktMTi2Pt)Ji{ri2kt) +PtJi{ri2Pt)MrM] H dacosT,2{a)kt 1, (7.12) 
T12 Jp^ kt J 

where the actual limits of integration over 77, 9, and a have yet to be put in agreement with the model we employ. 
Now, we have approached the most subtle point of our analysis. This expression includes triple integrations, any 
of which (if performed formally) yields singular functions. For the sake of definiteness, let us start with the second 



term in Eq.(7.12) (which corresponds to the forward scattering of soft quark on a hard gluon from the distribution 
ng{a, kt)), rewriting it in its most expanded form, 

■^a f 1 ai 2 N (■''1 ~ ■''2) cosh?7/2 



J2^—^ I dr]0{T^2) Jab'i2{v)Pt] / da J i[ti2 [r]) k t] cos[ktTi2 {r]) cosh{a - ^p{r]))] dh . (7.13) 

7ri?l J T12 J_oo Jp. 

The first observation is that at large kt (which is the condition that the distribution ng(a,kt) can be measured 
within a short time) the main contribution to the a-integration comes from the domain a « tp{r]) where the phase of 
the cosTi2{a)kt is stationary. This is the domain of collinear interaction when the hard gluon from the distribution 
Ugia, kt) has almost the same rapidity as the virtual quark with transverse momentum ql — kt + pt in the self-energy 
loop. Obviously, this quark is also hard. Furthermore, its propagator, G[,.et] (ti, r2; (/t) — 9{ti — T2)G'[o](ti, r2; gt), is 
devised only from the free on-mass-shell partial waves which themselves are well localized in the rapidity direction. 
Hence, we deal with the intuitively very clear case of collinear absorption and emission of the gauge field quantum 
by a charged particle. All participants of the process are moving with the same velocity. According to the property 
of localization of states in wedge dynamics studied in paper [II], such a fine tuning of a to V' is indeed possible. 
This is illustrated by the left-hand figure on FigJ|. where the grey segments of the hyperbolas t — T2 and r = n 
correspond to the rapidity intervals occupied by the soft quark mode, rpt < 1, at the beginning and at the end of the 
scattering process, respectively. The bold black and the dashed segments show the rapidity intervals where the hard 
virtual quark and the hard gluon are localized at the same times. All three fields effectively overlap around rapidity 
ri2 = — ?//2 at r = r2 and around rji — +J//2 at r = ri. The rapidity direction between these points is exactly ^{i]). 
The rapidity a of the external gluon is sufficiently small and is close to the rapidity ipirj). 



15 





FIG. 4. Geometry of fields in the forward scattering amplitude qg — > qg. 



The maximal rapidity width of the interaction domain is defined by the causality condition rj^2 > 0: which immedi- 
ately establishes the upper boundary |r/| < r^o- Since the coUinear interaction corresponds to the condition a w ^{ff}, 
the rapidity of the hard gluon must be within this geometrically defined interval as well. The opposite case is depicted 
in the right-hand figure of Fig.^. The rapidity tpi.'n) is so large, that the external gluon is not localized within the causal 
boundaries ±?7o/2 of the interaction domain. In order to avoid this, we have to impose an even stronger requirement 
that IV'I'?)! < '70- According to ( [7.8| ), we have |?A(?7)| > ??• Hence, we must further take \rj\ < 77*, where the boundary 
77* is defined by the equation ip{rj^) — 770, 

T1+T2 



which has a solution, 



tanh 



— tanh — 

T2 2 



tanh 7/0 



77* 



(7.14) 



(7.15) 



'1 -r '2 '''1'''2 

We remind the reader that ^ ^ 1; only this condition allows one to introduce the the time-dependent transverse mass 
/i(pt, r). In order to simplify further analysis, it is convenient to present the internal integral over kt as the difference, 



T12 



JiiTi2kt) cos[Ti2{a)kt] dkt , 



(7.16) 



where the first singular term is the integral over kt, computed from to 00, and thus, it completely accounts for the 
domain kt —f oo. It includes the function 



f{T1,a) = [Tl2{cc) 



'12J + 



^/^ = [rf2(77)sinh2(a- 7^(77))];^/^ , 



which is singular at a = 7/^(77), thus fully accounting for the expected coUinear enhancement. This function, however, 
is a canonical distribution with respect to both its arguments a and 77, and it carries the standard regulators for 
the subsequent integrations. We shall consider the singular and the regular terms separately. Using the above found 
limits, we may write the singular term as 



-'2 



, (n - T2) cosh 77/2 
drj- 

-»?. 



1 



n2(?7) 



Ti2{ri)\ 



da 



[sinh^ 



(7.17) 



16 



where, since Ti2Pt ^ 1, we put Jo[Ti2Pt] ~ 1- After an obvious change of variable, the internal integral of Eq. (7.17) 
can be split into two, 



Vo 



da 



-\a-ip\ 



-7)0 [sinh la - 
Since by the definition of the (+)-distribution. 



,1/2 



/3 



e"" da 



l3 



'da 



sinh a Jq sinh a 

we obtain the singular part in the form, 

'"'* (ti — T2) cosh 7^/2 drj 
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da 


.Jo 


'I 


sinh a 




da 




3-2/3 












Je " 
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jsing 



r In {- [1 + e-4w _2e-2'w cosh 2^/^(77)]} 



(7.18) 



(7.19) 



Next, it is convenient to trade 77 for a new variable y, Ti2{ri) = (ti — T2)y. The helpful relations for this change of 
variables are 

dy _2 , (ri - T2Y 



cos]i{r]/2)dri 



-1 



(ri-r2)2 



^12 1??*) 



1 + (ti - T2Y /tiT2 



2Tir2 



cosh 27/) = 



(T1+T2)" 1 
2Tir2 7/2 



2riT2 



Taking into account that cosh 27/; — > 1, when 7/ — > 1, we obtain, 



jSing 

-'2 — 



\AT^ Jy, yV^^y^ 



In- 



1 



-2?7o 



Tl - 7-2 Ti - T2 

In ■ 



(7.20) 



2 VT1T2 ^/TiT2 ^TxT2 

This formula has two distinctive elements. The first element is the large In^, which is due to the collinear geometry 
of the interaction. This would lead to a divergence if the interaction domain were unlimited. The second element 
is the small factor ^ which is due to the small volume occupied by the interaction and it completely suppresses the 
potential divergence. One may notice that when the mean time ^Jt\T2 increases, then ^ ^ 0, and the corresponding 
part of the self-energy also tends to zero. This is easy to understand, since with the mean time growing, the system 
becomes more and more diluted locally. 



The regular part of Eq. ( 7.13 ) is given by the integral. 



jreg 
^2 



dr] 



(n - T2) cosh 77/2 



da 



J\(jx2kt)co'&\Tx2{oi)'kt\ dkt , 



(7.21) 



where, when — t^q < 77 < t^q, then Ti2{a) varies between its minimal and maximal values. 

Therefore, when fcf < Pt, we have Ti2kt ~ {ti — T2)kt <^ 1, and both functions under the integral over fcj can be 
expanded in Taylor series. All integrations become trivial and yield R 



Even if we impose no limitations on a the estimate is still as small as 



n 



^2 



^/T^^2 



Tl - T2 



rj^ , ^7 22) 



We have chosen this form of the answer, because our major assumption is vaUd only as long as rpt < 1 and because 



Ja ^j:^ / dr] 2 Pt >/i[n2(?y)pt] 



the dispersion equation (B^) has a kinematic factor y^rvT^ in it. 

The other two terms in Eq. (7.12) can be studied along the same guidelines. The third term is suppressed with 
respect to J2 by the factor pt/h, which is small by our major model agreement and it could have been discarded on 
this ground only. To be on safe side, let us rewrite it as 

n- ^ (T1-T2) co sh 77/2 

1 - Jo[Ti2{-n)kt] cos{Ti2{a)kt) 

u,L, ^ - 7ii - ill L^^— e J -r / 
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where the first term is the integral over kt from to 00. As i could be expected, the integrand is regular. Since 
there is an accounted for difference between kt and qt^ the exactly collinear regime becomes impossible and we do not 
have the large collinear logarithm in J^. Overall, this term is also suppressed at least by a factor ^ stemming from 
Ji[Ti2(ri)pt\ in the integrand. 

The first term in Eq. ( 7.12|) corresponds to the forward quark-quark scattering with high momentum transfer, 

Ji = J dr] Ja[Ti2{'n)pt] 

^'^ ; 177 r , ,1 ^coshg / MTuiri)qt]sm{Ti2{e)qt) dqt \ . (7.24) 

Here, we again recognize the collinear singularity which is, as previously, regulated by the (+)-prcscription. All further 
calculations for jTi are similar to the case of J2 and the answer reads, 

js^ng^-S^ + ^ein^ ^ j-.^2^C^ (7.25) 

y/TlT2 VT1T2 

These results will serve for us as a reference point for the estimates of the mathematically more complicated part 
connected with the radiation field of the transverse magnetic TM-modes. Before we address this issue, it is expedient 
to look at the obtained results more attentively and trace the correspondence between the calculations and physical 
pictme in more details: 

1. It has been observed in Sec. ^ (for the vacuum part of the quark self-energy) that in the framework of wedge 
dynamics, the collinear problems do not jeopardize the field theory. In "material part" of the self-energy, the collinear 
interactions were proved to be the most intensive and to lead to a visible enhancement of the interaction between the 
quark and radiation field. However, this enhancement never turns into a disaster of collinear divergence. One of the 
trivial reasons is that the space-time domain of the interaction is now limited, and large logarithms are multiplied by 
small phase volumes. 

2. A deeper insight into the wedge dynamics, shows that even intermediate collinear singularities observed in 
the terms J7i and ^2 are spurious. In order to reveal this fact, one can notice that the singularity ai a — ipi^rj) is 
present only in J7i and 1/2- It is absent in J'3, because of the extra negative power of kt brought by the subleading 
term of the angular integral df. This extra k^^ effectively suppresses the distribution ng{a, kt) at large kt- Next, one 
may ask, what minimal change of rig (a, kt) at large kt is necessary in order that the intermediate collinear singularity 
does not appear at all. This can be learned by changing the order of integration in Eq. (7.13). One can start from 
the integral drj with an assumption that the integrand only slowly varies within some interval of a around a = 0. 
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Then, it is easy to see that the singular term J^™^ of Eq. (7.20) totally originates from the domain of kt oo. It 
comes from a logarithmic integral between two infinite limits. This residual piece emerges only because we extend 
the distribution Ug ~ k^^ (obtained from a dimensional estimate) to an arbitrarily large kt- As we have already 
mentioned, this dependence is unphysical, e.g., because the distribution cPkt/k^ is not normalizeable. It has to be 
modified above some value of kt and, therefore, the singular term must vanish completely. Thus, in wedge dynamics, 
the phenomenon of coUinear enhancement is intrinsically connected with the basic property of localization inherent in 
the one-particle states. Only the states with infinitely large kt can have a precisely given rapidity and be responsible 



for the singularities like we encounter in Eqs. (7.16) and (7.24). 

3. Our way to pick out the leading contributions (in the mixed representation of wedge dynamics) from the 
space-time domains, where the phases of the interacting fields are stationary, is a generalization of the known method 
of isolating the leading terms using the pinch-poles in the plane of complex energy. The similarity of two methods 
can be easily understood since, e.g., the quark density correlator in the self-energy can be presented as a sum of two 
propagators, G^i^{q) — G[oo](9) + G'[ii](g). In the plane of the complex energy — k^ + p°, the (Feynman-type) 
propagator G'[oo]('?) has poles in the second and fourth quadrants, while the (anti-Feynman-type) propagator G[ii]{q) 
has poles in the first and third quadrants. The radiation part of the retarded gauge field propagator D^ret] (k) has poles 
in the third and fourth quadrants. Therefore, in both terms of G[i]{q)D[j.et]{k) = G[oo] (9)£'[ret] (^) + G^n] (g)£'[ret] (fc), 
the integration path along the real axis of the complex k'^ plane is pinched between two poles (one from I?[rei] j ^^nd 
the second from G[oo] or Gjnj ) giving the leading contribution when p'^ is small, and the three-momenta k and 
q — k -\- p are are very close to each other. Similar arguments are valid for the second part, G[,.et]Z3[i] , of the quark 
self-energy. The term G[i] (5)-D[;o„g] (k) is exceptional, because the propagator -D[;ong] (^) of the longitudinal field has 
no poles corresponding to the propagation. 

The wedge dynamics does not allow for a standard momentum representation, since its geometric background is 
not homogeneous in the t- and z-directions; accordingly, we do not have familiar pinch-poles in our calculations. 
Nevertheless, the patches of phase space where the phases of certain field fragments are stationary and effectively 
overlap, do now the same job as the pinch-poles, and yield the same answers when the homogeneity required for 
the momentum representation is restored. The way wedge dynamics tackles the problem is genuinely more general, 
because it addresses the space-time picture of the interacting fields.^ The momentum space is now split into the 
subspaces of rapidity and transverse momentum; the correlation between the particle's rapidity and its location is 
increasing with the increase of its transverse momentum . The role of pinch-poles is taken over by the geometrical 
overlap of the field patterns with the same rapidity. This observation can serve as a footing for the future development 
of an effective technique for perturbative calculations in wedge dynamics. The arguments of the localization are not 
applicable to the longitudinal part of the gluon field. In the term G[i]I?[;o„g] , no patch in space-time is dynamically 
selected, since -D[;ong] is not assembled from the propagating waves that could match the virtual quark in the loop by 
their phase. This is in line with the absence of pinch-poles due to i'liong] {k) in the momentum picture. 

The arguments presented above allow one to estimate the contribution of the radiation fields of the TM-mode in 
a very economical way. Let us consider the group D^^g^^^ which is very similar to the term J7i studied above, as 



It is well known, that the threshold behavior of the imaginary part of the photon self-energy can be derived from the pinch 
geometry of the poles of the electron propagators Since at the threshold, the e^e~ pair is created with zero relative 

velocity, the pinch indeed corresponds to the overlap of the stationary phases of the e"*" and e~ wave functions in the maximal 
possible volume. 
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an example. Now, the invariant ^-'[q] , as can be inferred from Eq. (A4), is known only in the integral representation, 
and not in an analytic form. Let us, therefore, employ the analytic form of ^Pj^,, given by Eq. (A22). It is easy to see, 



that the integration over a in T^^^^ leads to the same causal step-function 0(t^2) and, as previously shown, we have 
1 77 1 < 770. The expanded form of this term is 



2 . (ti - T2) cosh 77/2 



T12 



da tanh (a — — ) tanh (a + — ^ 



Yi[Ti2{v)qt] sm[Ti2{a)qt]dqu (7.26) 



where we have integrated out the azimuthal angle (f assuming that kt,qt 3> pt (the first correction is smaller by the 
factor {pt/ktY <C 1). The internal integral over qt can be transformed into 



cosh(Q: — '0('7))] 



1/2 
+ 



Yi{Ti2qt) s\n[Ti2{a)kt\ dkt 



(7.27) 



which brings us very close to Eq. (7.24) for JJi. Once again, we encounter a collinear singularity at a = ip(r]), and 



exactly the same arguments force us to set the same limits in the integrals, as in Eq. ( 7.24 ). We do not have to 
continue the calculations to understand the smallness of J74, mention only, that due to the narrow limits of two 



rapidity integrations in Eq. ( 7.26 ), the product of the two hyperbolic tangents in the integrand will add extra to 
the order of smallness of J^4 . By the same argument as used previously, we can easily learn that the singular term in 



Eq. ( |7.27D is spurious. 

In this group, associated with the TAf-mode of the radiation field, the leading (still parametrically small, and equal 
to J7i) contribution comes from the D^''^ component of the gluon correlators. We may summarize by the observation 
that only the overlap of the domains of stationary phase in two correlators matters. It can be visualized via the 
partial-wave expansion of any of the correlators in the self-energy loop. 



VIII. NONLOCAL PART OF THE LONGITUDINAL PROPAGATOR IN THE MATERIAL PART OF 

THE SELF-ENERGY. 



The longitudinal part of the gluon propagator contributes, to the invariant E", the term 



S°(ri,r2) 



[long] 



An 



{ri-T2 



d^kt 



V 



(2) 
[Ion 



T1T2 



-V 



ivv) 

long] 



(8.1) 



Using Eqs. (AI6) and (A17) for the gluon invariants, and the first of Eqs. (A22) for the quark invariant g'^^^, we arrive 
at the following expression which accounts for the non-local part of the longitudinal propagator (the contact part was 
studied in Sec. 

[long] 



S"(ri,T2) 



sCfK 



271^ Rl 



^ 0(ti - T2) 



dqt 



di] 



(ti - T2)cosh77/2 
\r?2\'^' 



cosh?/ r\-tq,shMn\(i^ 



T1T2 



dt 



{tI^) Y,iT,2qt) + -0{-Ti2) K,if,2qt) 



(8.2) 



where we have integrated out the dependence on the azimuthal angle ip in the approximation of qt ^ pt- The key 
observation that allows one to judge about the smallness of this term is that the limits of integration over variable 
t are very close, and the factor [1 — t'^/TiT2] is very small. This factor reflects a known competition between the 
electric and magnetic interaction of moving charges which reduces the net yield almost to zero. Let us replace t by 
the dimensionless u = t/ ^Jt\T2- Then this part of the integral becomes 
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T1T2 



■ ■ - dt^ 



du 



Vl-?V4+C/2 



\/i-?V4-e/2 



(8.3) 



where the hmiting behavior /(^) const jS^ when ^ — > can be conjectured from the behavior of the functions ^1(2;) 
and K\[x) at small x. However, the Laplace transforms of these functions in Eq. ( |8.2| ) are singular functions of 77 and 
we have to be careful in estimating these terms. In fact, the ^ order of smallness is not altered by the remaining 
integrations. First, it is useful to notice that the last factor in square brackets in Eq. (B^) is nothing but the invariant 
(/j'j^j which, according to its integral representation (A22), equals zero at ^ = 0. Next, it is profitable to change the 
variables of integration in the following way. We trade 77 for y according to T\2{f]) = (^i — T2)y in the domain rj^2 > 0- 
In the complimentary domain = — ^12 < 0, we change 77 for y using ti2(?7) — (ti — T2)y. We also replace ^qt by a 
new variable q. 

iasCp-N'f 





[long] 








mat 




dy 






dy 





0{ti ~ T2) 



1 



qdq 



1 + f 



I Y.ir^qy) e'^-^^^^ ^KHi-y^)^ 



Ki{Tmqy) e 



(8.4) 



where, we remind the reader that Tj„ = \/t\T2. When the argument is small, the functions ^1(2;) and K\(x) are. 



2 X 



In 



X X 



2 V 2 



0{x^ 



Kiix) 



0{x^ 



- +ln- 

x 2 

It is easy to see now that in the sum of the two integrals dy over the interval (0, 1), the leading singularities dy/y 
exactly cancel each other. Furthermore, it is safe to take the limits of ^ — > 0, w — > 1 in the integrand, and even to set 
the lower limit of the integral dq to be zero. The resulting integral is convergent, ^-independent, and yields a term of 
the order In the remaining integral, the variable y runs from 1 to 00, the integrand is not singular at finite y and 
q, and it is exponentially suppressed at large y and q. The behavior of the integral at ^ — > is not singular and the 
basic upper estimate const x remains unchanged for the entire integral 



A). 



ACKNOWLEDGMENTS 

The authors are grateful to Berndt MuUer and Edward Shuryak for helpful discussions at various stages in the 
development of this work, and appreciate the help of Scott Payson who critically read the manuscript. 



IX. APPENDIX. WIGHTMAN FUNCTIONS AND PROPAGATORS OF WEDGE DYNAMICS. 

In this Appendix, we put all field correlators into a form which is needed for the practical calculation of the 
self-energy. The density of states Z?[i] and the causal part Dp^ of the gluon propagator are used in the form of 
decomposition over the transverse modes. 
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where 



rvy 





g-ifctT COSh(Q-?7). — 



1 



47r3/2fct 



-/2 



(A2) 



are the transverse electric and transverse magnetic modes of the radiation field found previously in paper [II]. Here, 
we denoted, 



fi{T,r)) =itanh(a-?7)(e-*'^'^™'''("-'') - 1) 



f2ir,v) 



-ikfT cosh(a—ri) 



+ iktT- 



cosh^(a-77) ' cosh(Q:-?7) 
Starting from this form, we get the components of the commutator _D[o] (t2, ti; 772 — 771; fct), ^ 



^[0] - 



da 
2^ 



^[0] 



kf kg 



da 



sin ktTi2 



k'p k g 



tanh (a + ^) tanh (a — ^) [ sin fct7i2 — sin ktTi + sin ktT2 ] , 
— [(1 + k^TiT2) sin ktTi2 - ktTi2 cos ktTi2 



(A3) 



(A4) 



2tt kf cosh" (a + f ) cosh" (a - f ) 

+ sin ktT2 — sin ktTi — ktT2 cos ktT2 + ktTi cos fcjTi ] , 



(A5) 



where Ti = n cosh(Q; — 77/2), r2 = T2 cosh(Q; + 77/2), ri2 = T1 — T2. In the first of these equations, the coefficients of the 
tensors {Srs — krkg/k'f) and krks/kf are the invariants and Pj^j^ of Eq. ( |2.9| ), respectively. The latter is due to 

the TM-mode of the radiation field. Up to the factor k^'^ , Eq. ( [A5| ) defines the invariant The underlined terms 

are connected with the boundary conditions imposed on the TAf-mode at r = 0. They cancel with the underlined 



terms in the longitudinal part of the gauge field propagator given by Eqs. (A14) and (A15) thus providing the causal 
behavior of the components -Dj^^^^jj and Dj^^^-^ of the retarded propagator -D[a(iii] (''2, ti). In the body of the paper, 
we call the residues of this cancelation as the diagonal components of D^. The "off-diagonal" components of the 



commutator ZJJqj are 



da tanh(a + ?) r 

7; TTT frr si'^ f'tTi2 - ktTi cos ktTi2 + sin fctTs - sin fc^Ti + ktTi cos ktTi 

27r cosh"(a - §) '- 



ikr f da tanh(Q; 



,,,, _ 
M ~ k? 



|- [ sin ktTi2 + ktT2 cos ktTi2 + sin ktT2 - sin ktTi - ktT2 cos ktT2] 



(A6) 
(A7) 



2tt cosh (a I 2 ; 

The commutator is not a symmetric tensor. However, by examination, these components are odd with respect to the 
rapidity difference ?y = 771 — 772, and hence they do not contribute to the effective quark mass we are computing in this 
paper. 

The tensor of the gluon density D^^, (t2 , ti ; 772 — Tii;kt) has the "diagonal" components. 





h h 







cos ktTi2 



k k T] 77 

- ^ „^ tanh(Q; H — ) tanh(a )(cosA:tTi2 — cosktTi — cosktT2 + 1) 

fcj 2 2 



(A8) 



'^In all formulae below, the gluon rapidity a is counted from the reference point (771 + ?72)/2, the geometric center of the 
coordinates rji and 772 of the vertices in the self-energy loop. Thus, it corresponds to the rapidity 6' in the integral representation 
of the quark correlators in paper [II] . 
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da 



W / 27r/c2cosh2(a + f)cosh2(a-f) 



— [(1 + fc^^iTa) cos ktTi2 + ktTi2 sin fctTia 



cos ktT2 — cos ktTi — ktT2 sin ktT2 — ktTi sin ktTi + l] 



(A9) 



from which one can infer the invariants and Pjj^j'' of Eq. (^^) exactly in the same way as it was done for the 

invariants of the commutator Z?[o] . The off-diagonal components, 

—kr f da tanh(a + |] 



l(2) 



1] kf J 27rcosh2(a- f) 



[ cos ktTi2 - ktTi sin fctTia - cos fcjTa - cos hTi - ktTi sin hTi + l] , (AlO) 



D 



.r,r _kr r da tenh^i 2XtcosktTi2 - ktT2smktTi2 - cosfctTs - cosfctTi - fctTssinfctTs + l] , (All) 



2tt cosh^(a + 2) 

are also non-symmetric and odd with respect to the rapidity difference rj. They also do not contribute to the effective 



quark mass. Equations (A8)- (All) give the components of the vacuum density of states of the gauge field in the 



wedge dynamics. In order to incorporate the "material" part given by the distribution of real gluons, the integrand 
of each of Eqs. (A8)- (All) must be multiplied by the common factor [1 + 2ng{kt, a)]. 



The full tensor of the longitudinal part of the propagator that defines the field A{ti) via the current j(r2) at all 
preceding times, 



1 [long] 



(ri) 



T2 dT2 d%Am"^V2, n; ?7l ^ kt)j"'{T2) 



was found in paper [III] in the following form, 
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dvd^k 

{2TT)^kl 



krks[Q-lA,y{kj_T2) - Q-l,t^{k±Ti)] krV[Qi,i^ikj_T2) - Q -l.i„{k±Ti)] 
vks[Q-l,i„{k±T2) - QlAi^ik^Ti)] l''^[Ql,iu{k±T2) - Ql,iu{k±Ti)] 



The diagonal components of this longitudinal part of the gluon propagator are just the differences between the vector 
potentials of the "static" gluon fields at the final time ri and the initial time T2, 



T^[long] _ krks \ COth |?/| , rtkt sinh |r)| _ -rafet sinh |7,K 
^rs 7,2 1 r, c I 



I — tanh(a + -) tanh(a - -) [ sin fc^Ti - sin fc^Ta ] I , (A14) 
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/cj sinh 1 77 1 kt sinh rj 2 sinh 1 77 1 
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2n kj cosh^(a + f ) cosh^(a 



— [ sin ktTi - sin /cfTa - ktTi cos ktTi + fcjTa cos /cfTa ] 



(A15) 



By the derivation, these components include 6(ti — T2) of the following origin. The source current which acts at the 
moment T2 produces the simultaneous longitudinal electric field E{t2). The gauge field potential is rebuilt from the 
electric field at the time ti > T2 by integrating the electric field E{t) over all times from to ti. The underlined terms 
in Eqs. ( A14 ) and ( [Al5| ) cancel out in the full assembly of the retarded propagator D[adv]{T2, ti) with the underlined 



terms in the radiation part, Eqs. (A4) and (A5). In the body of the paper, we call the residue of this cancelation as 
the diagonal components of which can be conveniently written as 
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j^[long\ _ _krks_ h COsh?] 1""^ ^-tkt sinh \ri\ 
H 2 J 



(A16) 



2 2 

Once again, the "off-diagonal" components of the longitudinal part of propagator, 
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(A17) 



[sin ktTi - sin fc^Ta + fc^Ta cos fctTa] } (A18) 
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signry 
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da tauh(n + ^) 
2^cosh"(a - f)' 



sin fctTi + sin fcjTa + fctTi cos fctTi] } (A19) 



are odd with respect to r] and do not contribute the effective quark mass. 

The gauge-field correlators have several distinctive features. First, the lengthy expression for each component is 
such that the gauge field correlators obey the boundary condition A^(r = 0, r*) = which provides continuity of the 
field at r = 0, and allows for a complete fixing of the gauge. Second, in the rs— and r?r?-components of the propagator 



D[Zv] ir2, TuV, h) = -e{ri - T2)D[^{T2,n,v; h) + £>f (r2, n, r?; h) 



(A20) 



the boundary terms cancel between the transverse and longitudinal parts. This fact provides causal behavior of the 
5](Ti,r2) that defines the dispersion law. 

The fermion invariants g]^a\ were derived in paper [II]. For the sake of completeness, we reproduce the final answers 
here, 
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